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Abstract. We investigate how to create entangled states of ultracold atoms trapped 
in optical lattices by dynamically manipulating the shape of the lattice potential. 
We consider an additional potential (the supcrlatticc) that allows both the splitting 
of each site into a double well potential, and the control of the height of potential 
barrier between sites. We use supcrlatticc manipulations to perform entangling 
operations between neighbouring qubits encoded on the Zeeman levels of the atoms 
without having to perform transfers between the different vibrational states of the 
atoms. Wc show how to use superlattices to engineer many-body entangled states 
resilient to collective dephasing noise. Also, we present a method to realize a 2D 
resource for measurement-based quantum computing via Bell-pair measurements. We 
analyze measurement networks that allow the execution of quantum algorithms while 
maintaining the resilience properties of the system throughout the computation. 
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1. Introduction 

The experimental realization of bosonic spinor condensates in optical lattices has opened 
up the possibility to exploit the spin degree of freedom for a wide range of applications 
[H [2]. Recently, further experiments have proved that these systems are promising 
candidates for quantum information processing purposes (see |3] and references therein) 
and the study of quantum magnetism [1]. In the literature there are several proposals 
using spinor condensates in optical lattice, e.g. to create macroscopic entangled states 
[5|, [6] or to explore magnetic quantum phases [7j. In this paper, we consider the 
possibility of using optical superlattices to manipulate the spin degrees of freedom of 
the atoms and engineer many-body entangled states. 

Superlattice setups allow the transformation of every site of the optical lattice 
into a double-well potential, and also the control of the potential barrier between 
sites [8], m H]. Atoms with overlapping motional wave functions interact via cold 
collisions that coherently modify their spin, while conserving the total magnetization 
of the system [TOl [TT] . Thus, by splitting each lattice site into a double- well potential, 
two atoms occupying the same site become separated by a potential barrier, and the 
interactions between them are switched off. Similarly, interactions between atoms can 
be switched on again for a certain time by lowering the potential barrier between the two 
sides of the double-well. Therefore, superlattice manipulations offer the possibility to 
apply operations between large numbers of atoms in parallel. In this work, we propose 
to exploit this parallelism to engineer highly-entangled many-body states. 

The process of splitting a Bose- Einstein condensate using double-well potentials has 
already been studied in the context of mean-field theory [121 [ISl [13 • ^^^^ P^^^ 

of this paper, we will derive a two-mode effective Hamiltonian that allows an accurate 
description of the system as the superlattice is progressively turned on. Our approach 
differs from previous treatments of double well potentials (see e.g. Ref. [I2]) or the work 
in Ref. [7], as the effective model we use is valid even for low barrier heights. We will use 
this effective Hamiltonian to analyze the process of splitting every site into a double- well 
potential starting with two atoms per site. We will show that this process creates a Bell 
state encoded on the Zeeman levels of the atoms in every site, each atom occupying one 
side of the double-well. Since these states are resilient to collective dephasing noise, we 
will use a lattice with a Bell pair in every site as a starting point to engineer many-body 
entangled states. 

In the second part, we present a new method for implementing an entangling 
\/SWAP gate that does not require to transfer the atoms between different vibrational 
states. We will show that the application of this entangling gate between (or inside) Bell 
pairs allows the creation of many-body entangled states resilient to collective dephasing 
noise. 

In the final part, we propose a new method involving superlattice manipulations 
to realize a 2D resource state for measurement-based quantum computating (MBQC) 
formally similar to a Bell-encoded cluster state. Since the realization of logical gates 
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between non-neighbouring qubits is practically very difficult in optical lattices [12], the 
one-way model for quantum computation is particularly relevant for these systems [161 
[TTlIlH]- The resource state we propose is resilient to collective dephasing noise, which 
makes it less prone to decoherence than the usual cluster states. Its utilization as a 
resource for MBQC requires adjustments of the measurement patterns used for cluster 
states, which we describe in the last section. We note that the realization of Bell- 
encoded cluster states as a resource for MBQC was proposed in Ref. [19] using lattice 
manipulations in three dimensions. The distinguishing feature of our proposal is that the 
resource state is created via 2D superlattice operations which have been demonstrated 
in the lab already. 

This paper is organised as follows. In Sec. [2] we derive an effective model for 
describing the dynamics of atoms within one site in the presence of a superlattice 
potential. In Sec. [3l starting with a lattice where every site contains two atoms of 
opposite spin, we show how to create a Bell pair in each site by splitting it using 
a superlattice potential. In Sec. HI we present how to perform a gate between two 
neighbouring qubits by manipulating the superlattice potential, and in Sec. [5] and Sec. [6] 
we show how to use this gate to create and probe many-body entangled states. Finally, in 
Sec. [7]we propose a method for creating a resource for MBQC via Bell-pair measurements 
and provide measurements network to implement a universal set of gates. We conclude 



We consider a gas of ultracold alkali atoms loaded into a deep optical lattice. The atoms 
are assumed to be bosons initialized in the / = 1 hyperfine manifold [2]. The lattice 
potential is given by VbL(r, s, 6) = VT[sm'^{nz/a) + sin^^ny/a) + Vs{x, s, 9)] where Vr is 
the depth of the potential and 



The depth of the potential will be given in units of the recoil energy Er = 
h^Ti"^ / (2Ma^) where M is the mass of the atoms. The parameter s G [0, 1] is determined 
by the relative intensities of the two pairs of lasers. The constant a corresponds to the 
size of a unit cell when s = (see Fig. [1^). We will refer to unit cells when s = as 
to lattice sites. Changing the value of s from to 1 transforms each site in a double 
well potential. We will refer to each side of this double well potential when s = 1 as to 
a subsite. The angle 6 allows for the manipulation of the potential barrier inside each 
site {6 = 0) or between sites {9 = 71/2). A few lattice profiles corresponding to different 
parameters s and 9 are shown in Fig. [H In the following, we will refer to the lattice 
profiles corresponding to the values of s ~ and s ~ 1 as to the large and small lattice 
limits, respectively. 

We assume the lattice depth Vt to be sufficiently large so that hopping can be 
neglected in the small and in the large lattice limits [20] . 



in Sec. E 



2. Model 




(1) 
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Figure 1. (a) The system is initialized with two atoms per site in the state 
Wa(r)^ (8) S (l/i = 1, TOi = +1; /2 = 1,7712 — —1)) and the lattice parameters are s = 
and 6^0. (b) The lattice potential is smoothly altered by increasing the lattice 
parameter s while 9 = 0. (c) At the end of the first step, the lattice parameter is 

5 = 1. Each lattice site then contains a Bell pair, (d-f) After setting the angle to 

6 = 7r/2, the other set of barriers arc lowered by decreasing the lattice parameter s 
from s = 1 (d) to < s < 1 (e) and back to s = 1 (f). 



The Hamiltonian of the system is given by 

H = Hk + Hz + Hint (2) 

where Hk describes the kinetic energy, Hz is the Zeeman term, and Hmt describes the 
interactions between particles. These terms are defined by [211 122] 

Hk = Uv J2 ^i{r)hoiv,s,e)^^{v), (3) 



-1,0,1 



Hz=[dr ^EzAB) ^i{r)^Ar), (4) 

o-=-l,0,l 

Hint = -] dr [co4o(r)ioo(r) + Cs iL(r)^m(r)J , (5) 



m=-2 

where cq = 47r/i^(2a2 + ao)/(3M) and C2 = 47rfi^(a2 — ao)/(3M) with ap the s- 
wave scattering length of the channel associated with the total angular momentum F, 
ho = {—h?'V^/2M) + VoL{'^, s, 6) and AEz^a{B) is the energy shift caused by the Zeeman 
effect in the presence of a magnetic field B = (0, 0, B) oriented in the ^-direction. The 
operator Apmpi'^) is given by 

/ 

Apmpir) = (F,mi.|/i,mi;/2,m2)^„i(r)^^^2(r), (6) 

mi,m2=-f 

where |/i, mi; /2, m2) = |/i,mi) |/2,m2) and (F, mjT'l/i, mi; /2, mg) is a Clebsch- 
Gordan coefficient. 
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In Appendix A, we derive the effective Hamiltonian describing the dynamics of the 
atoms at one site using the field operator 

¥,{r)=aiwa{r) + blwbir), (7) 

where aj-(fej-) is the operator that creates a particle with spin \ f = l,mf = a) {a = 
{— 1, 0, 1} denotes the Zeeman level) in a motional state associated with the symmetric 
(anti-symmetric) mode function Wa{wh)- The mode functions are centered in the middle 
of the site, and their shape depends on the lattice parameter s. The field operators 
obey the canonical bosonic commutation relations [^o-(r), ^)^/(r')] = 5crcr'i5(r — r') and 
[^o-(r), "^a'iT^')] = [^t(r), ^)^/(r')] = 0. We will see in the next section that the inclusion 
of a second motional mode in the Hamiltonian allows an accurate description of the 
system's dynamics in both the large and the small lattice limit [23]. A more detailed 
explanation of the interaction term ([5]), as well as the full Hamiltonian of the system in 



terms of creation (annihilation) operators can be found in Appendix A 



3. Creation of a Bell state on every lattice site 

In this section, we use the effective Hamiltonian introduced in the previous section 
to show how to generate a lattice where every site contains a Bell state encoded on 
the Zeeman levels of the atoms. This procedure will be used as a preliminary step to 
engineer several types of many-body entangled states. 

We start with a deep optical lattice in the large lattice limit (LLL) and two atoms 
per lattice site in the state 

IV^init) = Wa{rY ® 5 (l/i = 1, mi = +1; /2 = 1, m2 = -1)) , (8) 

where S is the symmetrization operator. In general, the two atoms can undergo spin- 
changing collisions coupling to various hyperfine levels. As we want to implement qubits 
on the Zeeman levels, we will limit the spin dynamics effectively to two hyperfine states. 
There are several ways to achieve this in practice: (i) One possibility is to exploit the 
conservation of angular momentum and choose two states such as |/ = 2, mj = +2) and 
1/ = IjTUf = +1), which are not coupled to any other hyperfine states by the interatomic 
interaction, (ii) An accidental degeneracy in ^''Rb offers yet another possible route: since 
the two s-wave scattering lengths oq and 02 are almost equal, transitions between the 
state 1/1 = l,mi = 0; /2 = l,m2 = 0) and states where particles have opposite spins 
occur on a time scale (cq + C2)/c2 ~ 802/(02 — cto) ~ 300 times slower than any other 
allowed transition. Hence, after initializing the system at time t = in state (IHl), we 
can limit our dynamical description to transitions between states with opposite spins, 
as long as the manipulations we are about to propose happen on a faster time scale 
than this transition. Notice that the use of fast-switched microwave fields that suppress 
spin-changing collisions [1] allows the relaxation of the constraint to be faster than the 
original time-scale of spin-changing collisions. The preparation of the state ([8]) has 
already been experimentally achieved [TT| [2]. 
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The first step in our proposal consists of raising tlie superlattice potential, 
i.e. changing the lattice parameter s{t) from the LLL at time t = to the small lattice 
limit (SLL) at some time t = T sufficiently slowly so that the ramp does not create 
excitations in the system. Once the superlattice is fully ramped up, each site is split 
in two, atoms no longer interact with each other and their spin state remains frozen in 
time. 

The shape of the mode functions Wa{r) and ^^^(r) depends on the lattice parameter. 
In the LLL, they correspond to the ground and first excited state of each lattice site. 
However, as the lattice parameter s approaches the SLL, the mode functions transform 
into the symmetric and anti-symmetric superpositions ii'a(r) = + v5/?(r)]/ -\/2 and 

Wb{r) = [v5L(r) — V^_R(r)]/v^, where the mode functions v^L(r) and </?R(r) are centered on 
either side of the double-well potential [23l[T2]. Hence, when the mode functions WL,R{r) 
become centered within each subsite [23], it is convenient to write the Hamiltonian of 
the system in terms of the operators 



li-bi 



(9) 



which create a particle with spin a on the left and right side of the double- well potential, 
respectively. In the basis 



Ti,) 
T,i) 
i,T) 

,n) = 



vac) 
vac) 
vac) 



(a{ + 6|)(a{ + fe{)/2|vac). 



{a\ + b\){d\ - 6{)/2|vac), 



6|)(a| + 6l)/2|vac), 



d\ -b\)id\ -6{)/2|vac), 



(10) 



the Hamiltonian of the system reads (see Appendix A ) 



H 



n 



( E + tJn 




Jn 


Xjl ^ 


Jji 


E + xn 


^ Xn 


Jji 


Jn 


Xjl 


E + xn 


Jn 


\ Xn 


Jn 


Jn 


E + Uu J 



where ^ = K + H with K = / dr <(r)/jo w^{r) and Kq = {-h'^V^/2M) + VoL{r, s, 6) 
is the sum of the single-particle energies, U^i = 'Jn^Uaa + Ubb + 6f/afe)/4 with 'j'\iUuu' = 
7lj / dr (|wi,(r)||u;i,/(r)|)^ and ■Jn = Cq — C2 is the generalized on-site interaction, 
Jn ~ (ya — Vb) f^+lui^aa — Ubb) /4: the generalized tunneling matrix element between the 
two sides of the well, and xn ~ 7TJ.(^aa+f^bb— 2f/a6)/4 corresponds to the density-density 
interaction energy between the two sides of the double-well. The latter is proportional 
to the overlap between the functions v^L(r) and (fiR^r). It cancels in the small lattice 
limit as a consequence of the localization of the mode functions v^l(i') and ^PR{r) inside 
each subsite. 

Our two-mode model does not only describe the system well in both limits, it also 
takes into account the effects of density-density interactions between two subsites when 
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Figure 2. Numerical calculation of the anti- fidelity 1 — !F where = |(?/>|?/>o)P with 
lip) the state of the two atoms at the end of the ramping. We have used the parameters 
of a lattice with period a = A/2, A = 840 nm, and Vt = 60Er, initialized with two 
^''Rb atoms per site in the state ([8]) at time T = and a magnetic field of -B = 30 G. 



the potential barrier between them is still not fully ramped up. This property makes it 
particularly valuable for the study of the system dynamics between the two limits. 
In the SLL, the Hamiltonian (fTT]) reduces (up to the energy shift E) to 



fysLL 



( Un 


-J 


-J 
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-J 








-J 


-J 








-J 
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-J 


-J 





(12) 



where Va = E — J and Vb = E + J with E = J drip}^{r)hoipL{r) is the single-particle 
energy and J = f dr (p}^{r)ho(pR{r) the hopping integral, and Ui;i = 7|; J dr \(Pl\'^ is the 
on-site interaction |23j . 

For zero tunneling J = 0, the Hamiltonian H^f'^ has two degenerate ground states: 
I 1, 1) and I i, t)- This degenercy is lifted at finite J where the symmetric superposition 



|V;o) = ^(IT,i) + U,T)) 



(13) 



is the ground and the antisymmetric superposition the first excited state. These two 
low-lying energy levels with an energy splitting of 4 J^/f/^j are separated from the other 
excited states by the interaction energy f/||. 

We have carried out a dynamical simulation of the splitting dynamics using the 
exact full Hamiltonian defined in Eq. ([2]) and Appendix A for two particles and one 
site. For the time-scales considered throughout this paper, non-adiabatic effects due to 
the changes of shape of the mode functions can be safely neglected [23]. We have used 
s{t) = sin^[t/(2T)], and computed the anti-fidelity 1 — JF with JF = |(^/;|^/;o)P between 
the state of the system at the end of the numerical simulation of the ramp and 
the Bell state ( IT3!) . For the parameters considered, we have found that the anti- fidelity 
reaches 1 — JF ^ 10~^ for a quench time of T ~ 1.1 ± 0.1 ms (see Fig. [2]). It is expected 
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that better fidelity can be obtained for larger values of U-n [23] . This observation makes 
our scheme relevant for current experimental setups. 

Thus, changing the topology of the optical lattice from the LLL to the SLL drives 
the state of the atoms within each lattice site from ([8]) to the maximally entangled Bell 
state (fT3l) . Non-adiabatic transitions to excited states are suppressed by opposite parity 
in the case of the first excited state and by the on-site interaction f/|| in the case of the 
other excited states. After this first step has been completed, a system of initially N 
sites contains k = 2N subsites and its state is given by a tensor product of Bell pairs 

N 

l<^>f ) = (g) l^o). (14) 

i=l 

Remarkably, the system in state (fT4l) is resilient to dephasing noise for a magnetic 
field slowly varying over a distance of one lattice site [3]. 



4. Implementation of an entangling gate 

Since the interactions between two atoms depend on the overlap between their wave 
functions, they can be dynamically switched on and off by lowering and raising the 
potential barrier between the two subsites. This is done by varying the lattice parameter 
s from s = 1 at time t = Otos = l— ?7 and back to s = 1 at time t = t. In the SLL, 
where J-n/Un <^ 1, tunneling can be treated perturbatively, and the Hamiltonian (fTTl) 
can be projected onto the subspace of singly-occupied sites. 

We distinguish two cases: either the neighbouring subsites are occupied by atoms 
with opposite or equal spins. For two atoms with opposite spins, we can use Eq. (ITT!) 
and the effective Hamiltonian in the basis | t, J.) and | J,, t) reads (up to a constant 
energy shift) in second-order perturbation theory 

»r=(f ..,.)(;;) 

Since the effective Hamiltonian commutes with itself at different times, we obtain the 
following analytical expression for the time evolution operator 

J:^(cff)(J^^ ( e-'-^cos^ -ie-'-^sin^ \ 



+ Xn]dt. (17) 




9 P 



u 



While the first term in Eq. (fTTl) is due to second-order tunneling [71 H] , the second term 
is due to density-density interactions. This term is negligible in the SLL, but when 
the potential barrier between the two sides of the double well is reduced, it adds an 
extra phase between neighbouring particles with different spins because of the difference 
between the scattering lengths Cq and C2- 
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For atoms with equal spins a = {t? 1} we start by writing down the basis states as 



I'^f^) = 4,<x4,<x|vac) 
in which the Hamiltonian reads 



vac) 

(4 + &l)(4-&l)/2|vac) 
{al-bl){dl-bl)/{2V2)\y^c), 



\ 



J. 

Xo 



Jrr 



Jrr 



Xo 



Xo 

J. 



\ 



J 



(19) 



where jaa = cq + C2 is the interaction constant for particles with equal spin a, 
Efja = K + Vfe ± Ez the sum of single-particle energies, Ez = —gjiBBr /2h the 
Zeeman energy shift in the first order, Jo-o- = {Va — Vb)/\/2 + ^aoiUaa — Ubb)/^ the 
generalized tunneling matrix element, Uaa = laaiUaa + Ubb + QUab) /'^ the generalized on- 
site interaction and Xo-o- = loaiUaa + Ubb — ^Uab) / ^ the off-site density-density interaction. 
In the SLL this reduces (up to the energy shift of -Eo-o-) to 

/ U^^ -V2J \ 
-V2J 
f/„ 



H. 



SLL 



\ 



(20) 



For the state |cr, a) there are no other low-energy states accessable, so that it acquires 
only the phase factor e"^ during the sweep. In second-order perturbation theory we find 
that K is similar to ( flTl) 



1 



4J. 



2 



+ Xo 



dt. 



(21) 



Since - J, -\/2J and '^^hoo = (co - C2)/(co + C2) ^ 1, we get k ^ 2</). 

Putting the results of this section together we find that lowering and raising the 
potential barrier between two subsites n and n + 1 implements the two-qubit gate U 



U 



( |oo) \ 















\ 


/ |oo) \ 


|01) 







e"^ cos (f) 


— ie"^ sin </> 







|01) 


|10) 







— ie^"^ sin 


e'"^ cos 







|10) 




v 













/ 


V 111) / 



(22) 



where \ij) = \i)n ® \j)n+i with |1)„ = cl^^\va.c) and |0)n = cj^^Jvac) where c|^_^ and c]^^-^ 
are the creation operators of a particle at subsite n with spin | and t, respectively. 

Numerical simulations of the system dynamics have been carried out using the exact 
full Hamiltonian for two particles and one site — that is, including the contribution of spin 
changing collisions to the nif = state — for a system of '^^Rb atoms, Vr = QOEji and 
rj = 0.3. We have computed the overlap between the states resulting from the numerical 
simulation with their analytical approximation and found that for the parameters 



Robust entangled states and MBQC using optical superlattices 



10 




Figure 3. Schmidt measure of the state ([29| as a function of T and k. The Schmidt 
decomposition of the state (|29l) was calculated by partitioning the system in equal 
parts, which yields the maximum possible Schmidt rank. 



considered and switching times r up to tens of milliseconds, the time evolution operator 
of the system is accurately approximated by Eq. (|22l) . For longer times r, the evolution 
of the system in the reduced Hilbert becomes non-unitary due to slow transitions to the 
state where the two particles have spins a = 0, as expected. 

Since the application of the gate U is realized in parallel on all qubits in the x- 
direction, the phase due to the Zeeman shift in the first order in B cancels for systems 
with an equal number of atoms in the a = +1 and a = —1 state. For such systems, we 
find that for ramping times r such that 0(r) = {7r/4, 7r/2, 37r/4}, lattice manipulations 
realize the gates (in the usual computational basis) [25] 

t/| = VSWAP (23) 
= SWAP, (24) 

Ua. = SWAP VS WAP = VSWAP^ (25) 

In the remainder of this paper we will show that, using a lattice of Bell pairs as a 
starting point, this gat^ can be used to create both many-body entangled states and a 
resource state for MBQC. 



5. Creation and detection of a state with a tunable amount of entanglement 

Depending on whether the lattice parameter ^ is set to ^ = or ^ = 7r/2, the gate 
Uff, in Eq. (122|) is applied between neighbouring subsites inside or between lattice sites, 

I In the context of optical lattices, different methods to implement this gate have been put forward for 
different encodings of the logical qubits [261 [27l [28] . Recently, a vSWAP gate between qubits encoded 
on the internal spin state of the atoms was realized experimentally by manipulating both the spin and 
motional degrees of freedom of the atoms by means of optical superlattices [3j . 
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respectively. Up to an irrelevant global phase, these operations read 

N 

l^^'ifi) = (g) t> (26) 

n=l 

or 

N-l 

W^"(f/) = i®(g) i. (27) 

n=l 

Using this notation, we define a knitting operator by 

t{U)=W\U)U'''"{U). (28) 

Starting from a lattice of Bell pairs, the state resulting from F successive applications 
of the operator /C(W|) on the state dHj) is denoted by 

\^f^)=mif\^f)- (29) 

In Fig. [3], we have plotted the Schmidt measure of the state fl29|) for up to 
= 20 qubits and partitions of size L = 10 as a function of F [291 [16]. We find 
that the Schmidt rank of the state is directly proportional to F, i.e. applying the 
entangling knitting operator (l28ll first connects the initial Bell pairs and then further 
increases the entanglement content of the state. The largest possible Schmidt measure 
raa.yi[P^{\^f^))\ = k/2 is reached after F = N/2 applications of the operator fl28|) . 
Once the maximum value for the Schmidt measure is reached, it remains unaffected by 
further applications of the operator (1281) . We note that since the operator U conserves 
the total number of spins, the state (l29l) is resilient to collective dephasing noise. Recent 
experimental results suggest that resilience to this type of noise significantly increases 
the decoherence time of the system [3]. 

Since the application of the gate /C(f/|) affects the density correlations in the lattice, 
its effect on the system can be observed experimentally via state-selective measurement 
of the quasi-momentum distribution (QMD) of atoms with spin up [301 E] 

2N 

" 2N 

As an illustration, we have plotted in Fig. |4] the QMD of a system of A; = 14 atoms 
after F = 1 and F = 7 applications of the operator /C(t/|) on the state fll4l) . Hence, 
starting from a lattice of Bell pairs, superlattice manipulations allow the creation of an 
entangled state with a tunable Schmidt rank that has a distinct experimental signature 
and is resilient to collective dephasing noise. These properties make this state suitable 
for experimental studies of many-body entanglement. 

6. Creation of maximally entangled states 

Together with site-selective single-qubit operations applied in parallel on every second 
atom, superlattice manipulations can be used to implement the entangling phase-gate 
operation 

Ci = U^{Z 1)Ue{1® Z), (31) 



^ E ^"^^'^^"^■^^''(4t^..t)- (30) 
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Figure 4. (a) Quasi-momcntum distribution of fc = 14 atoms after the application 
of the gate IC{Ujl)^ on the state dUD for (a) T = 1 and (b) T = 7. After T = 7 
apphcations of the knitting operator, the state is maximaUy entangled. 

where Ci = diag(l, — i, — i, 1) and Z is the Pauh matrix defined as Z = diag(l,— 1) 
in the usual computational basis. The operation Ci is related to the operation 
Cz = diag(l, 1, 1, — 1) by applying the single-qubit gate ® on every site of 
the lattice. A proposal for realizing arbitrary single-qubit gates on individual atoms 
in an optical lattice can be found in Ref. [31]. Notice that the single-qubit operations 
required to implement the gate (l3T|l between sites can be performed in parallel, and so 
each site need not be addressed separately. The implementation of a phase gate directly 
followed by a SWAP gate 

(SWAF) Ci = Ui{Z (g) 1)U^{1 ® Z) (32) 

is accomplished by adjusting the switching time r of the last gate. The operation ( l32ll 
is equivalent to (SWAP) Cz up to the unitary transformation ® \/^. Since the 
Cz operations between different qubits commute, k/2 successive applications of the 
operator K:{{SWAP)Cz) on the state (|±) = (|0) ± \ l))/V2) create a complete 

graph state \Kn). Complete graphs have the property that each vertex is connected 
to all the vertices of the graph, i.e. a complete graph of k vertices contains k{k — l)/2 
edges. Graph states of complete graphs are equivalent to the maximally entangled state 
\GHZ) = (10)®'^ + |l)®'=)/\/2 up to local unitary operations. 

Superlattice manipulations alone allow the realization of the state 

|$|-Ph) =^bw(f;t)^^^t)(fc/2)-l|^(0)^ (33) 

4 4 

resulting from k — 1 successive applications of the operator ( l25|) between and inside each 
site. The process leading to state ( l33l) is schematically represented in Fig. [5l This state 
is resilient to collective dephasing noise, and possesses a structure similar to a complete 
graph state (see Fig. [5]). Via the brute- force numerical calculation of the parameters 
of single-qubit unitary operations, we have found that fl33|) is locally equivalent to a 
complete graph state of the same size for up to k = {4, 6, 8} qubits. This suggests that 
this property holds for an arbitrary number of qubits. 
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J 






Figure 5. Schematic representations of tlie process leading to state (|33|) . Qubits 
initially forming a Bell pair are linked by a dark blue line. The boxes depict the 
location where the VsWAP^ = SWAPVSWAP is performed. After a VSWAP 
gate has been applied between two qubits, they are connected by a (light) blue line. 
Numerical calculations show that for up to 8 qubits, the state ((55)) is equivalent to a 
complete graph state under local unitary operations. 

The resilience of this state to collective dephasing noise, and its symmetry 
properties make it a good candidate for the improvement of the sensitivity of quantum 
spectroscopic measurements in noisy environments (see Ref. |32j). 



7. Creation of a resource for MBQC 

In this section we will show how to create a state useful for MBQC via the application 
of the Cz gate via lattice manipulations in both the x and y direction. This state is 
formally similar to a Bell-encoded cluster state, and hence its utilization as a resource for 
MBQC only requires an adjustment of the measurement networks used for cluster states. 
To demonstrate the universality of our resource for MBQC and derive the measurement 
networks required to perform quantum algorithms, we will employ the method recently 
developed by Gross and Eisert in Ref. [33l [M], which connects the matrix product 
representation (MPR) of a state with its computational power. A review of the basics 



concepts presented in Ref. [33] can be found in [Appendix B 

The MPR for a chain of k systems of dimension d = 2 (qubits) is given by 

I^a)= E {R\At^■■■At^\L)\H■■■^,). (34) 

ji...ifc={0,l} 

It is specified by a set of 2k D x D-matrices, which we will refer to as the correlation 
matrices, and two D-dimensional vectors \L) and \R) representing boundary conditions. 
The parameter D is proportional to the amount of correlation between two consecutive 
blocks of the chain. Notice that the right boundary condition vector \R) appears on the 
left. This choice improves the clarity of calculations later when we will use the graphical 



notation explained in Appendix B.2 
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Starting from a lattice of Bell pairs, we find that the state resulting from the 

(0)\ 

k I 
. ^[«fe-ii 



application of the gate lA^^ 



{Cz) on the state \^\ 



n...ifc={0,i} 



H]) has the MPS 



(35) 



where 



m 



|+)(0|, iw = |-)(i|, 

ii)(oi, m =|o)(i|, 

\R) =v^|0). 

Here, each atom is labelled from 1 to /c, and the correlation matrices A'^/^J and 

are associated with odd and even atoms, respectively. Equivalently, it can be written as 



(36) 
(37) 
(38) 



\^c) 



1 



(i^|C<[n]...C<M|^)|-^...-^)^ 



(39) 



where 





ii...ijv=0 






= it^iira = 


|-)(o| 




= iM^ii] = 




|0) 


= |io), |i) 


= |oi) 



(40) 



In this representation, each site is labelled from 1 to and the correlation matrix C^^/^^ 
is associated with a pair of atoms. The measurement of odd and even atoms in the X- 
eigenbasis Bx = {|±) = (|0) ± |l))/-\/2} on the state (135|) projects the auxiliary matrix 
associated with the measured atom onto the eigenstate (|0) + (— 1)*|1))/ \/2 depending on 
the measurement outcome s (see Eq. flB.2p in Appendix B] ). In the graphical notation, 
this implements the operations 

(41) 



A[X] oc X'H, 



B[X] 



oc XZ' 



(42) 



for odd and even atoms, respectively. Here, the operator H is the Hadamard gate 
and X the Pauli-X matrix [25]. Similarly, measurements in the 0-eigenbasis = 
{(|0)±e^<^|l))/v^} yield 

- ocX'HS{(j)), (43) 
► ocXZ'S{(f)), (44) 



B 



where >S'(<^) = diag(l,e''^) in the usual computational basis. Also, we find that since 
the measurement of the zth and {i + l)th atoms {i odd) in the Bell basis has only two 
possible outcomes, we have 

(45) 



C[X] 



oc ZX'H, 
oc ZX'HS{(j)), 



(46) 
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A D ,4 D 




X 



Figure 6. Starting witli ID system (cc-direction) in tlie state ([35|l . a 2D state useful 
for MBQC is obtained by applying Cz gates between even subsites (denoted by the 
letter B) in the y-direction. 



where the measurement eigenbases are given by Bx = {(|0) ± \l))/\/2} and = 
{(|0) ± e''^|l))/v/2}, respectively. Since any single-qubit gate can be decomposed into 
Euler angles, i.e. 



U,o,{C,v,0 = S{OHS{v)HS{0, 
the application of one of the following sequence of measurements 









Am 








A[X,] 


-^B[4>,]^A[Xe]^ 
















C[X,] 

















(47) 

(48) 
(49) 



on the state (13511 realize arbitrary qubit gates, up to a known by-product operator Uy.. 
For instance, applying the measurement sequence fHHj) with measurement outcomes 
s = (1, 0, 0, 0, 1, 1), where Si denotes the outcome of Xi or (pi, implements the operation 
f^sf^rot(02, 03, 05) with f/j] = HZX ou the state of the correlation space. Similarly, 
the measurement sequence ( 149!) with outcomes s = (0, 0, 1, 0) implements the operation 
?^sfA-ot(02, 03, 04) with f/s = ZXZ. Although both measurement sequences (HHl) and 
( H9|) preserve the system's immunity to collective dephasing, less measurements are 
required using the sequence ( l48l) . The state of the correlation space can be read-out 



using a scheme described in Appendix B 



Superlattice manipulations can also be used to engineer 2D systems. In the 
remainder of this section, we will consider the MPR of the 2D state \ip2D) resulting 
from the coupling of ID systems in the state (l35l) via the application of a Cz gate 
between even subsite^ in the y-direction (see Fig. [6]). The MPR of this state is given 
by 



-)r(Ob® |+)„(0|d. 



(50) 



§ See e.g. Refs. [35l |3] for relevant 2D superlattice setups 
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IT 



\+)r{l\l(S)\-)u{l\d- 



(51) 



Also, the expansion coefficients of |^/'2d) are represented by 



(7[«2,1 

I 



R 



D 



I 



R 



(52) 



D 



wliere \U) = |0) and |D) = |+). 

This state is formally similar to a cluster state where qubits are encoded on two 
atoms and the |+) state corresponds to the Bell state f|T3|) . Although the application 
of the Cz gate between sites in the ^/-direction implements the same operation between 
encoded qubits, it corresponds to a different operation between encoded qubits in the x- 
direction. Therefore, it is not possible to execute algorithms designed for cluster states 
using the state \tp2D) as a resource by simply interchanging, e.g. X measurements with 
Bell-pair measurements X. 

In the remainder of this section, we will show that \ip2D) constitute a universal 
resource for MBQC by presenting how to control information flows through the lattice, 
and by providing a measurement network that implements an entangling two-qubit gate. 

Since the tensors of Eqs. fl50|) and flSTl) factor, they can be represented graphically 

by 



Using this representation and Eqs. (!50|) and (1511) . we find that 



(53) 



C[Zi] 



C[X2] 



oc XZ'^+'^+'^'H. 



(54) 



Thus, measurements in the basis Bx cause the information to flow from left to right, 
and measurements of the vertically adjacent sites in the Bz basis shields the information 
from the rest of the lattice [33]. Since we can isolate horizontal lines, they can be used 
as logical qubits. 

An entangling two-qubit gate between two horizontal lines is realized as follows. 
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— C[X5] 

, ~cz , 



(55) 



c)^C[X,] 



where the middle site is measured in the basis By = {(|0) ± i|l))/v^} and c G {0, 1}. 
The lower part of the network reduces to 



Jl 



c)-C[Xi 



oc 



± 



Z'\+) 



(56) 



-iy^xH\c) 



Plugging (!56|) into the middle line yields 



C[Z^]\^C[Y2]\^C[Z^] 



oc 



(57) 



Finally, plugging (IFH) into the upper of the network gives 



oc Uyi-iZ) 



(58) 



where f/^ = e'5(iX)**H^XS'(|)^Z'*5(— iZ)''2+'^3+i_ Thus, the measurement network flS^ 
implements an entangling controlled-phase gate between the upper and lower horizontal 
lines up to a know by-product operator. Since arbitrary single-qubit gates can be applied 
on each horizontal lines independently, this completes the proof of universality. The way 
of dealing with by-product operators at the end of the measurement sequences does not 
differ from the usual MBQC scheme with cluster states (see e.g. [IE])- Notice that since 
single and two-qubit gates are performed via the application of pairwise measurements, 
the system remains invariant to collective dephasing at any time during the execution 
of an algorithm. 



8. Summary 

We have analyzed the dynamics of a bosonic spinor condensate in a superlattice 
potential, and shown that a lattice with a Bell pair in every lattice site can be realized 
via the dynamical splitting of each site. We have proposed a scheme that allows 
the application of an entangling vSWAP gate between and inside Bell pairs. The 
successive application of this gate between and inside lattice sites was shown to create 
an entangled state with a tunable Schmidt rank that is resilient to collective dephasing 
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noise; and a maximally entangled state which, as numerical evidence suggests, is locally 
equivalent to a GHZ state. Finally, we have presented a state that is obtained by 
connecting Bell pairs in two dimensions via an entangling phase gate, and shown that 
it constitutes a resource for MBQC formally similar to a Bell-encoded cluster state. We 
have provided measurement networks for implementing a two-qubit entangling gate as 
well as arbitrary local unitary operations. Our implementation has the advantage that 
it allows the execution of quantum algorithms while leaving the system unaffected by 
collective dephasing noise. 
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Appendix 

Appendix A. The full Hamiltonian 

Using the field operator (I7j), the single-particle terms of the Hamiltonian ([2]) in second 
quantized form are given by 

-^0 = Hk + Hz 

= [Va + AEzAB)]aia^+ [Vb + AEzAbMK, (A.l) 

o-=-l,0,l cr=-l,0,l 

where = J dTwl{r)hoWu{r). Assuming that the magnetic field is weak {B < 200 G), 
the Zeeman shift AEz^a{B) is accurately approximated to the second order in B [36] 



if a = -1, 

AEzAB) = I -^-^1^ ifcr = 0, (A.2) 

^bB 3{gfiBBf .„ 
r-^; 11 (T = i, 

4 leA^hf 

where g is the gyromagnetic factor, Ai^hf is the hyperfine splitting energy, and /^^ is 
the Bohr magneton. 

In the low-energy limit, atomic interactions are well approximated by s-wave 
scattering, with the scattering length depending on the spin state of the colliding atoms. 
At ultracold temperatures, two colliding atoms in the lower hyperfine level /low will 
remain in the same multiplet, since the interaction process is not energetic enough 
to promote either atom to a higher hyperfine level /high [in]- Also, since alkalis have 
only two hyperfine multiplets, the conservation of the total angular momentum by the 
scattering process implies that interatomic interactions conserve the hyperfine spin /i 
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and /2 of the individual atoms. Thus, the interaction between two particles located at 
positions ri and r2 is given by p!Ol 



Mnt(ri - ra) 



M 



-6{ri - rs) J2 ^P^f 

F=0,2 



(A.3) 



where F = /i + /2 is the total angular momentum of the scattered pair, is the 
projection operator for the total angular momentum F, and ap is the s-wave scattering 
length of the channel associated with the total angular momentum F. The operator 
is given by [37] 

F 



^ \F,mF){F,mF\ 



(A.4) 

mp=—F 

where |F, m^?) is the state formed by two atoms with total angular momentum F and 
niF = rrii + m2 with mi and m2 the projection on the quantization axis of /i and /2, 
respectively. The quantization axis is defined by the direction of the magnetic field 
present in the system. Boson statistics require the state |F, mi?) to be invariant under 
particle exchange, and hence only the terms corresponding to values oiF = {0, 2} appear 
in the sum of Eq. (1A.3|) . Using the relation fi ■ f2 = ^2 — 2^0) where fj = (/*, /*, /*) is 
the total angular momentum operator of the particle i, the interaction operator can be 
re-written as [ini [38] 

(A.5) 



Mnt(ri - r2) = 5{vi- Y2) [cq + ■ f; 

From Eq. (lA.Sj) we see that interactions conserve the quantum number niF [36| [39]. 
Using the field operator ([70, the interaction term in second quantized form reads 
1 



int 




Co 



dridr2 6{ri - r2) x 



*(r2 

*t(r2)/2*(r2; 



+ C2 

£=x,y,z 

where : 6 : represents the operator o in normal ordering, \E'(r 
and the matrices fg are given by [22] 

/ 1 \ 

1 1 
1 



1 

72 



\ 







ri 
J V 



I 



V2i 



1 







1 
1 
-1 





/ 


1 





\ 


' J z 













I 


\ 








-1/ 



(A.6) 

(vl>„i(r),vl>o(r),^i(r)) 
(A.7) 



The Hamiltonian given in ([5]) is a compact form of (]A.6p . In its expanded form, Eq. (]A.6P 
reads 



-2 J"' 

(Co + C2) f^I(r)^l(r)^i(r)^i(r) + ^li(r)^li(r)^_i(r)^„i(r) 
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+ 2(co + C2) [^I(r)^/,(r)^i(r)^o(r) + ^U(r)^|,(r)^_i(r)^o(r) 
+ 2c2 [^S(r)^S(r)^i(r)^_i(r) + ^l(r)^Li(r)^o(r)^o(r) 
+ C0 [^S(r)^|,(r)^o(r)^o(r) 

+ 2(co - C2) [^I(r)^li(r)^i(r)^„i(r)J |. (A.8) 

The interaction Hamiltonian (1A.8I) can be expressed in terms of cij. and 6j. operators 
using the relation 



(^a;^aj^6^/6y + hlh\a„ia^^ + (^aj.6jy + hlal^ (a^'by + j ,(A.9) 
where U^,y' = J dr {\w^{r)\\w^i{r)\)'^ and a, 7 = -1, 0, +1. 

Appendix B. Measurement-based quantum computations and MPSs 

Appendix B.l. General idea 

In this section, we present the general idea developed in Ref. [31]. From Eq. flM|) . one 
can see that if the ik site in the computational basis and the outcome is obtained, 
then the state \L) of the auxiliary system becomes 

|L)' = A^1|L). (B.l) 

Thus, measurements change the state \L) of the auxiliary system. Similarly, if the 
measurement of a local observable at site k yields an outcome corresponding to the 
observable's eigenvector then the matrix Aj!*^ transforms into 

A[0], = (0,|O)a[°1 + (0,|1)aW, (B.2) 

and the vector of the auxiliary system becomes 

\L)' = m,\L). (B.3) 

From this point of view, a measurement on some physical site change the correlation 
properties between this site and the rest of the chain, i.e. it performs an operation on 
the auxiliary system state |L), which we will sometimes refer to as the state of the 
correlation space. 

Appendix B.2. Graphical notation 

In order to deal with the MPS of higher- dimensional states, the graphical notation 
introduced in Ref. [M] is very helpful. In this notation, tensors are represented by 
boxes, and their indices by edges. Vectors and matrices are thus symbolized by 

1^)=^, (B.4) 
{R\ = , (B.5) 

(B.6) 



aI^ . 
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Consequently, vector and matrix operations can be represented graphically by 
{R\L) - 



R 



BA 
BA\L) 



A^B 



L^A 



B 



(B.7) 
(B.8) 
(B.9) 



Appendix B.3. Read-out scheme 



In order to use the state of the correlation space to process quantum information, we 
must be able to read it out at the end of the computation. It turns out that in our 
case a measurement of the {i — l)th site in the computational basis corresponds to a 
measurement of the correlation system just after the ii\i site. To illustrate this fact, 
assume that just after the measurement of the zth site, the correlation system is in the 
state |0), that is 

■ ■ ■ ■ ■ ■ - = |0), 



B 



(B.IO) 



or 



L 



B[<P,]^A[<P2] 



Thus, using Eqs. (l36l) and (1371) we have 



0) 




A[l] 


— oc|-)(l|0)=0. 


|0) 




B\l] 


^ oc |0)(1|0) = 0. 



=io). (B.ii) 

(B.12) 
(B.13) 

Hence, the probability of obtaining the result 1 for a measurement on the (i — l)th site 
is zero. Consequently, if the state of the correlation system is in state |0) after the ith 
site, then the {i — l)th physical site must also be in that state. Since 

oc |+)(1|0) = 0, (B.14) 

the same observation applies if one measures the {i — 2)th and (z — l)th sites [i odd) 
in the basis Ez = {|0), |1)}. Therefore, as a similar argument applies if the state of the 
correlation system is in the state |1) after the ith site, the description of the read-out 
scheme is complete. 



|0)-C[1] 
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